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Abstract. Bremsstrahlung of W and Z gauge bosons, or photons, can be an important 
dark matter annihilation channel. In many popular models in which the annihilation to 
a pair of light fermions is helicity suppressed, these bremsstrahlung processes can lift the 
suppression and thus become the dominant annihilation channels. The resulting dark matter 
annihilation products contain a large, energetic, neutrino component. We consider solar 
WIMP annihilation in the case where electroweak bremsstrahlung dominates, and calculate 
the resulting neutrino spectra. The flux consists of primary neutrinos produced in processes 
such as XX ~~ ^ vvZ and XX ~ > v&W, and secondary neutrinos produced via the decays of 
gauge bosons and charged leptons. After dealing with the neutrino propagation and flavour 
evolution in the Sun, we consider the prospects for detection in neutrino experiments on 
Earth. We compare our signal with that for annihilation to W + W~, and show that, for a 
given annihilation rate, the bremsstrahlung annihilation channel produces a larger signal by 
a factor of a few. 
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1 Introduction 

High-energy neutrinos produced by dark matter (DM) self-annihilation in the Sun offer an 
appealing signal for indirect DM detection [1-8]. See, e.g., refs. [9-15] for recent work on this 
topic. DM particles which pass through the Sun can scatter, lose energy, and become captured 
by the gravitational field of the Sun. They accumulate in the Sun's core and eventually reach 
a density where self-annihilation becomes important. Most Standard Model (SM) particles 
produced in the annihilations interact in the Sun and are hence absorbed, while neutrinos 
will escape to be detected in experiments on Earth. The neutrinos may be produced directly 
in the DM annihilations, or as secondaries in the decays of other SM particles. A high energy 
neutrino signal from the Sun would thus be a compelling indication of dark matter. 

The strength of this signal is controlled by the DM-nucleon scattering cross section 
which controls the DM capture rate and, if in equilibrium, the annihilation rate. Existing 
limits have been placed by Super-Kamiokande [16], IceCube [17-19], and other experiments. 
The strength of these limits depends on the annihilation channel assumed for dark matter, 
with limits for soft annihilation channels such as \X ~^ bb being weaker than those for hard 
channels such as XX ~^ W + W~ . A comprehensive study of the neutrino spectra for the 
main annihilation channels with 2-body final states was performed in ref. [9]. In the present 
article we study DM annihilation processes with 3-body final states, namely bremsstrahlung 
processes in which a W or Z [20, 21] or 7 [22-27] is radiated. In certain models, these 
bremsstrahlung processes have been shown to be the dominant DM annihilation mode. 

For non-relativistic dark matter, it is useful to parametrise the annihilation cross section 

as 

a AV = a + bv 2 + 0{v A ), (1.1) 
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where v is the dark matter relative velocity (in units of c). The a term arises only from 
s-wave annihilation, and the L th partial wave contribution is suppressed by a factor v 2L . For 
highly non-relativistic DM with wCl, s-wave annihilations, if permitted, are thus expected 
to dominate. In the galactic halo, v 2 ~ 10 -6 ; v is even smaller for the DM accumulated 
in the solar core. However, in many Majorana DM models, s-wave annihilation to a pair 
of fermions, XX ~^ f f ■, ls helicity suppressed by a factor of (mf/M x ) 2 due to a mismatch 
between the net helicity and chirality of the fermion pair [28, 29]. A well known example is 
the annihilation of neutralino dark matter in supersymmetric models. 

The circumstances in which annihilation to a pair of fermions is velocity and/or helicity 
suppressed were extensively detailed in ref. [21]. It is well known that radiation of a photon 
is able to lift this suppression via the electromagnetic (EM) bremsstrahlung process XX ~ > 
//7 [22-27]. However, as shown in refs. [20, 21], electroweak (EW) bremsstrahlung is also 
able to lift this helicity suppression, such that the XX ~^ ff^, XX ~^ fvW, and XX ~^ f fl 
annihilation channels all dominate over XX ~^ ff- Physically, the suppression is lifted because 
the emitted gauge boson carries away a unit of angular momentum, allowing a fermion spin- 
flip, so there is no longer a mismatch between the chirality of the leptons and their allowed 
2-particle spin state. Relative to the 2-body process, the rate for the 3-body final state suffers 
a suppression factor from the radiating particle given by a/ir. However, this may be more 
than compensated for by the consequent removal of the (rrif/M x ) 2 suppression factor. 

The importance of these radiative corrections to the dark matter annihilation process 
have been examined in a number of recent papers [30-41], with a focus on the indirect 
detection of dark matter due to annihilations in the Galactic halo. Interestingly, a possible 
hint of an EM bremsstrahlung signal was found in recent Fermi data [42]. In the present 
paper we consider the indirect detection of dark matter via annihilation in the Sun, noting 
that the EW bremsstrahlung annihilation channels produce large, energetic, neutrino fluxes. 

This paper aims to produce the final neutrino energy spectra, as they would be detected 
on Earth, for solar WIMP annihilation via the EW bremsstrahlung process, and to compare 
the resulting signal with that for other commonly assumed annihilation channels. In section 
2 we describe the neutrino spectra at production, including primary neutrinos produced 
directly in the annihilations, and secondaries produced from the decay of emitted gauge 
bosons. In section 3 we consider the scattering and flavour evolution of the neutrinos as they 
propagate through the Sun, to determine the final spectra, while in section 4 we compare 
the event rates for bremsstrahlung annihilation with those for other DM annihilation modes. 
Our results are summarised briefly in section 5. 

2 Neutrino production 

During the process of capture in the Sun, DM particles scatter elastically with solar nuclei, 
losing enough energy to become gravitationally bound. The capture rate, C®, is proportional 
to the DM/solar nucleon scattering cross section and depends upon astrophysical parameters 
such as the local DM velocity and density. Once captured, subsequent scattering causes the 
DM to become concentrated in the solar core, in a region of size ~ 0.01-R© ^100 GeV/M x , 
where Rq is the solar radius; the region is sufficiently small that we may assume all DM 
annihilations occur in the centre of the Sun [9]. 
The differential neutrino flux is given by 




(2-1) 



A: 
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where d is the distance from the centre of the Sun to the detector, and BR^ is the branching 
ratio to a specific final state, k. The annihilation rate r ann depends on the capture rate C® 
as 

r ann = ^C tanh 2 (t o /T4), (2.2) 

where to = 4.5 Gyr is the age of the Sun and ta is the time-scale for competing capture 
and annihilation processes [6, 9, 43-45]. This time scale is given by t a = 1/VC®A®, where 
^4® = (av) /Ves is the annihilation cross section times the relative WIMP velocity per volume 
(see, e.g. the reviews in [46, 47] for explicit expressions for these rates). This time-scale is 
often much smaller than to, such that an equilibrium is reached between the capture and 
annihilation processes with r ann ~ C®/2. In this case, it is the DM-nucleon scattering cross 
section, rather than the annihilation cross section, which determines the annihilation rate. 
We shall return to the question of capture- annihilation equilibrium in section 4.2. 

2.1 Dark matter coupling 

We use the leptophilic Majorana DM coupling of ref. [48] as a simple example of a scenario in 
which EW bremsstrahlung dominates the annihilation. The DM is a gauge-singlet Majorana 
fermion which annihilates to leptons via the interaction term 

f(u r) L e K\ x + h.c. = f(u L r, -l L r, + )x + h.c. , (2.3) 

where / is a coupling constant, e is the 2x2 SU(2) invariant antisymmetric matrix, and 
(?7 + , 77 ) form a new SU(2) doublet scalar which mediates the annihilation. An identical 
coupling arises in supersymmetric models if we take x to be a (bino-like) neutralino and 77 a 
sfermion doublet. 

Ref. [48] proposed this model as an explanation of the sharp excess in the cosmic ray 
e + /(e~ + e + ) fraction at energies beyond ~10 GeV, as recently measured by the PAMELA 
experiment [49-51]. Here we adopt this model as a simple prototype of the situation in which 
bremsstrahlung processes dominate the DM annihilation. As explained in [20, 21], the s-wave 
contributions to the lowest order annihilation processes \X ~~ vv, e + e~ are suppressed by fac- 
tors of (mi/M x ) 2 . Bremsstrahlung lifts this suppression to become the dominant annihilation 
process. 

We shall consider three possibilities for flavour structure of the branching ratios: DM 
annihilates to all three lepton flavours equally (/ e : : f T = 1 : 1 : 1), to electron-flavoured 
leptons only (/ e : : f T = 1 : : 0), or to /i and r flavours equally (/ e : : f T = : 1 : 1). 
However, we shall see that neutrino mixing reduces the sensitivity of the final results to the 
flavour structure of the couplings. 

2.2 Neutrino spectra at production 

We determine the energy spectra for the neutrinos at production by adding the contributions 
from the primary neutrinos produced directly in DM annihilation, and the secondary neutrino 
produced by decays of gauge bosons or charged leptons. Throughout, we neglect neutrinos 
with energies less than ~ 0.5 GeV. Note that IceCube, for example, has a detection threshold 
of ~ 100 GeV [52], which is reduced to ~ 10 GeV with the addition of DeepCore [53]. We shall 
be particularly interested in the higher-energy regions of the spectra, as the bremsstrahlung 
process leads to neutrino spectra with a distinct high energy peak; moreover, detection cross 
sections grow (approximately linearly) with energy. 
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Figure 1. The primary spectra of muon-flavoured neutrinos for /j, = M^/M^ — 1.01 (dashed, blue), 
/i = 1.2 (solid, purple), fi = 2 (red, dotted), /i = 5 (orange, dot-dashed) and fi = 10 (green, small- 
dashed). Increasing /i reduces the high-energy peak but increases the number of mid- range energy 
neutrinos. We have taken the annihilations to have equal branching ratios to all lepton flavours, 
f e = ffx = f T - Note that the x-axis is defined as x = E v jM x . 



2.2.1 Primary neutrino spectra 

We first consider the primary neutrinos, those produced from the direct DM annihilation 
channels: 

XX ->■ l~viW + , 

XX -> l + ihW~ , (2-4) 
xx -»■ fJ^-Z • 

The spectra of primary neutrinos from these processes are given in [20] , along with calcula- 
tions for the relative cross sections. The contribution from each bremsstrahlung channel is 
included according to these relative cross sections; for example the branching ratio for emis- 
sion of a W~ is given by a ev w /ctotaij where utotai is found by adding all the bremsstrahlung 
components together (including the contributions from XX ~ > Z and XX ~> ^~7> which 
do not produce primary neutrinos). Annihilation directly to neutrinos via the suppressed 
process XX ~* would produce a monoenergetic spectrum with E v = M x ; we neglect the 
contribution of this subdominant process in our analysis. 

The primary neutrino spectrum is sensitive to the ratio of the DM and propagator 

M 2 

masses, /i = We focus on small values of \x as this maximises the branching ratio of the 

X 

bremsstrahlung processes relative to the 2-body final state annihilation processes; however, 
bremsstrahlung remains the dominant annihilation channel for \i < 10. In figure 1 we plot 
the primary neutrino spectra for various choices of fj,. The spectra have a broad energy 
distribution, and a distinctive high energy peak at x ~ 1. Increasing /i decreases the height 
of this peak, but increases the number of neutrinos of mid-range energy. For the remainder 
of this paper we shall assume a value of \x = 1.2 (but will also show results for \x =5 in table 

1)- 

In figure 2(a) the primary spectra of muon neutrinos are indicated by the dashed green 
line. Here we assume equal branching to all lepton flavours, and so the primary spectra are 
flavour- independent; changing the branching ratios will clearly affect the flavour dependence. 
Note that the primary neutrinos are the dominant contribution to the total neutrino spectrum 
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Figure 2. a) The contributions to the injected spectra, assuming branching to all lepton flavours 
equally, from primary neutrinos (green, dashed), secondary neutrinos from gauge boson decay (purple, 
dot-dashed) and secondary neutrinos from decay of primary taus (orange, dotted); the solid blue line 
shows the total, b) The total injected spectra for different branching ratios: annihilation to all lepton 
flavours equally (black, solid), annihilation to e-navoured leptons only (blue, dashed), and annihilation 
to /i and r flavours equally (red, dotted). All anti-neutrino spectra are identical to the corresponding 
neutrino spectra, and so are not shown. Note that the s-axis is defined as x = E u /M x , and that we 
are using /i = 1.2. 

for each flavour (shown as the solid blue line in figure 2(a)), particularly for the high-energy 
region of interest. Comparing the spectra for different choices of the DM mass, we see that 
the primary spectra maintain the same general shape. However, increasing the DM mass 
relative to My//z decreases the effect of kinematic thresholds and moves the peak closer to 
x = 1. The antineutrino spectra are identical to the corresponding neutrino spectra. 

2.2.2 Secondary neutrino spectra 

We next consider the secondary neutrinos that are produced from gauge boson decays, or 
the decays of primary tau leptons. The decay of the emitted gauge boson is modelled using 
the Monte Carlo simulator PYTHIA [54], in which we create a W or Z at rest which is then 
allowed to decay with constraints to be outlined below. The neutrino spectra from the gauge 
boson decays are then boosted into the lab frame as per the appendix of ref. [21]. In the 
Sun, the core density is p = 140 g cm -3 , so energy loss processes are relatively important 
as the DM annihilation products may interact with solar material and lose energy before 
decaying further. Adopting a procedure similar to that in ref. [9], the relevant energy losses 
are accounted for as described below. 

Gauge bosons: Once produced, the lifetimes of the W and Z bosons are short enough 
that they decay before interacting and losing any significant amount of energy. Therefore, in 
PYTHIA, we set these bosons to decay as if in vacuum. The products of the hadronic decay 
modes of the gauge bosons are dealt with as explained below. 

Charged leptons: The typical stopping time of a charged lepton is T stop ~ 10~ n s. If this 
is shorter than the boosted lifetime, r stop < 77d ec , the lepton will lose all its energy before 
decaying. For the energy range of interest, E < 1 TeV, the muons are stopped before they 
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decay while taus decay without losing a significant amount of energy. Thus, in the PYTHIA 
routine, we discard electrons and muons, and allow taus to decay as if in vacuum. 

Quarks: Heavy hadrons have comparable r stop and 7Tdeo of order 10~ n s and 10 _12 s 
respectively, so will often lose some energy before decaying; light hadrons are usually stopped 
before decay. As a reasonable, simple, approximation, we allow heavy hadrons (those includ- 
ing charm, top and bottom quarks) to decay, while light hadrons are discarded. 

The resulting neutrino spectra from and Z decay, as seen in the lab frame, are 
shown by the dot-dashed purple line in figure 2(a). These secondary neutrinos pile up and 
become dominant at low energies, as they are often produced through a chain of decays 
starting at the gauge boson. For this case of even branching ratios (/ e = = f T ) the 
secondary spectra of v e and are equal, as light hadrons and muons are stopped before 
decay, while r, c and b produce equal numbers of each. As with the primary neutrinos, the 
antineutrino spectra are identical to their corresponding neutrino spectra. 

In the case where the coupling f T is non-zero, there will be a primary tau component 
which will decay to contribute to the secondary neutrino flux. As with the gauge bosons 
decays, the secondary neutrino spectra from tau decay (again created using PYTHIA) must 
be boosted into the lab frame and convolved with the primary tau spectrum. The neutrino 
spectra resulting from this process are shown by the dotted orange line in figure 2(a). The 
tau decays make a minor contribution to the v e and i/„ fluxes and a larger, though still 
subdominant, contribution to the v T flux. 

Adding the contributions of primary and secondary neutrino fluxes, we plot the total 
spectra at production as the solid blue line in figure 2(a). The total spectra are also shown in 
figure 2(b) for varied branching ratios. We note that the primary neutrinos make the domi- 
nant contribution to the total neutrino flux (summed over all flavours), while for particular 
choices of the flavour branching ratios the secondary neutrinos can dominate the flux for an 
individual flavour. However, once neutrino mixing is taken into account the primary neutri- 
nos will make the dominant contribution to the flux for each flavour. Therefore, our results 
are robust against uncertainties in the calculation of the secondary neutrinos, and instead 
are determined mainly by the flux of primary neutrinos (whose spectra can be calculated 
analytically) . 

3 Propagation in the Sun 

Following production, the neutrino flux is modified by flavour oscillations, and by scattering 
with solar nuclei, which results in some absorption of the flux at high energy. These effects are 
well-understood, and various techniques exist to deal with them; we shall follow the method 
of ref. [55] to account for the scattering interaction, and then adopt results of ref. [56] when 
evaluating the flavour evolution. Many authors combine scattering and mixing effects into 
a single 3x3 density matrix equation that describes the neutrino evolution. However in the 
energy range of interest, E v < 1 TeV, oscillations and scattering approximately decouple [55], 
as the matter effects in the Sun ensure negligible mixing between v e and i/^ T until after 
scattering becomes unimportant. We assume throughout that the v„ and v T fluxes are equal, 
as Un and v T mix maximally (to good approximation) and, for the energy range of interest, 
the oscillation length is less than or approximately equal to the interaction length [55]. 
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3.1 Scattering and absorption 

The evolution equation, used to evolve the neutrino flux from production to r c (the radius 
where interactions become unimportant), has the form 



dpi _ dpi 
dr dr 



dpi 
■ ■ dr 

inj 



| dpi 

NC dr 



(3.1) 

cc 



for neutrinos of flavour I, where p{r,E)dE is the neutrino flux. The first term on the right 
corresponds to the neutrino injection due to the annihilation of DM (including all secondary 
neutrinos). The second and third terms correspond to neutrino scattering through neutral 
current and charged current interactions respectively. 

The probability of an interaction between a neutrino and a solar nucleon is given by 



int 







R& a tot (E)N s (r')dr' = |, (3.2) 



where citot is the neutrino total interaction cross section and Ns(r) ~ Noe~ r / KR ® is the 
nucleon density of the Sun, with k = 0.1 and Nq = 1.3 x 10 26 cm~ 3 . The constants £ 
parametrise the energy scale at which interactions become significant, and are given by 

S v ^ 140 GeV, £ v ^ 213 GeV. (3.3) 

Note that the total interaction cross section combines neutral current (NC) scatterings 
v\N -f-7- v\N and i>iN o v\N (effectively removing a neutrino from the flux and re-injecting 
it at lower energy), and charged current (CC) scatterings v\N l~N' and ViN l + N' 
(effectively removing a neutrino from the flux and producing an almost collinear charged 
lepton). If I = e, p, then the final lepton loses most, if not all, of its energy before it can 
decay again, and so the neutrino is considered absorbed. If I = r, the resulting decays 
promptly, either leptonically or hadronically. 

Defining the optical depth as 

y ^ = W^lT f N s( r ')dr', (3.4) 



the evolution equation (3.1) takes the form 



ay h{E) + e 



poo 

E Pj (y,E)+ dE' 9j (E/E>) Pj (y, E) 
Je 



(3.5) 



The term h(E) represents our injected flux, the second term accounts for absorption of the 
flux through scattering, and the last term the reinjection or regeneration of neutrinos at 
lower energy. At high energies, absorption through scattering is a significant effect. Tau 
regeneration is also important, as absorbed neutrinos are repopulated and pile-up at lower 
energies. Again, this effect is most significant for large DM mass. The regeneration functions 
<7i are given, e.g., in refs. [9, 55]. Following ref. [55] we use a perturbative solution of eq. (3.5) 
to handle the regeneration. For simplicity, we take gi = 0.5 and have checked that for 
the largest M x we consider (where regeneration is expected to occur most strongly) this 
simplification produces an error <20% in the final spectra for energies > 100 GeV, which 
does not qualitatively change our conclusions. 
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3.2 Flavour oscillations 

We now include the effect of neutrino flavour mixing. We use the neutrino mixing parameters 
given in ref. [57], and assume the normal mass hierarchy, where mi < m,2 < m^. We shall 
adopt a non-zero value for #13, as indicated by recents results from the T2K [58] and the 
Daya Bay [59] collaborations, and demonstrate how the results change as we vary the value 
of this parameter. 

Neutrinos undergo flavour conversion during their propagation through the Sun. These 
flavour effects differ from oscillations in vacuum due to the presence of neutrino refractive 
indices. In the dense solar medium, the electron neutrinos experience a matter potential of 
V e (r) = v / 2G'pA r e (r), where iV e is the electron number density. At high density, this leads to 
an approximate alignment between the flavour and effective mass eigenstates, and suppresses 
oscillations of v e with f^, r - The matter potential decreases as the neutrinos propagate out- 
ward from the centre of the Sun, and the neutrinos undergo an MSW resonance [60, 61] 
when 

Am 2 /2E ~ V e . (3.6) 

At low energies, the propagation through this resonance is adiabatic and the neutrinos remain 
in given mass eigenstates. At higher energies the evolution becomes non- adiabatic, with 
the neutrinos experiencing a level-crossing, or jump between mass eigenstates. Assuming 
an exponentially varying solar profile (p(r) oc e _r / r °), with ro = -Rq/10.54, the crossing 
probability is conventionally given by [62] 

g 7Cos 2 6» _ j 

Pc = e7 _ x , 7 = 4vrr A , (3.7) 

where A = Am 2 /AE U . For antineutrinos, the v e potential is —V e {r), which is equivalent to 
making the replacement 6 — > | — 9 [63], thus 

P c = P c (9^^-9) . (3.8) 

We shall be concerned with two level crossings, controlled by Am^ and A.m 2 2 respec- 
tively. As Am| 2 is almost two orders of magnitude larger than Am^ the crossings decouple 
and can be treated independently. This separation is useful as it allows us to use a simple 
two-neutrino formalism which we apply to each crossing independently. This leads to 




P= P l c 1-P^0 01-P£ P£ , (3.9) 

V 1/ \o P h c l-P&J 

where Pjk represents the probability of a neutrino that originated in matter mass eigenstate 
k ending up in vacuum mass eigenstate j. Here, P\~, and Pq correspond to the lower and 
higher level crossing probabilities, driven by #12 and #13 respectively. See ref. [64] for a 
detailed discussion. 

We shall need to evaluate the level crossing probabilities for the two resonances, and 
for both neutrinos and antineutrinos. In some of these cases the resonance condition (3.6) is 
never actually met; nonetheless, the level crossing probability is non-zero and the expression 
given m eq. (3.7) still holds [56, 65, 66]. 1 For some parameters of interest (namely, high 



As discussed in refs. [56, 66], the crossing probability is largest when adiabaticity is maximally violated, 
and is non-zero even away from resonance. 
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Figure 3. The level-crossing probabilities Pc\ our chosen value of 6*13 = 8.8° is used in fig (a), 
while fig (b) shows the corresponding result for 6*13 = 0.1°. The first (^-driven) level crossing 
probabilities (Pq) are shown by the red dot-dashed (neutrino) and dotted (antineutrino) lines, while 
the second (#i2-driven) crossing probabilities (Pq) are shown by the blue bold (neutrino) and dashed 
(antineutrino) lines. 



neutrino energy) the lower level crossing takes place in the convective zone of the Sun where 
the solar density profile deviates from the exponential form assumed in eq. (3.7). We account 
for this by following the procedure developed in ref. [56], where Pq is calculated according 
to eq. (3.7), but with tq replaced by an energy-dependent numerical fit. 

The level crossing probabilities, for neutrinos and antineutrinos, are plotted as a function 
of neutrino energy in figure 3. Figure 3(a) shows the probabilities used here, assuming 
#13 = 8.8°, as suggested by recent T2K [58] and the Daya Bay [59] results, while figure 3(b) 
shows the probabilities if we instead allow On = 0.1°. Note that for very small #13, Pq for 
neutrinos approaches 1 for all energies (i.e. the resonance is completely non-adiabatic). This 
is to be expected because as #13 — > the heaviest vacuum eigenstate becomes completely 
decoupled, and does not mix with the other mass eigenstates. For antineutrinos Pq ~ 0, as 
the resonance condition is never met. We also see that the lower crossing probabilities, P l c 
are comparatively larger for neutrinos than antineutrinos, as v e is produced in the lightest 
matter mass eigenstate and, again, does not undergo a resonance. 

We apply eq. (3.9) to obtain the neutrino fluxes at the edge of the Sun, in the vacuum 
mass basis. The neutrinos propagate from the Sun to the Earth in these vacuum states, and 
can then be converted back to the flavour basis. The final spectra are shown in figure 4. 

3.3 Final spectra 

The final neutrino spectra, per annihilation, for the electroweak bremsstrahlung annihilation 
mode are shown in figure 4, for various choices of the branching ratios to the lepton flavours. 
Also shown, for comparison, are the spectra that would result from DM annihiation to 
W + W~ and t + t~ final states. These comparison spectra were calculated using the same 
procedure described in section 2. 

We note that while the initial spectra at production were sensitive to flavour branching 
ratios chosen, neutrino flavour mixing in the Sun considerably reduces these differences. 
For example, for branching only to the electron flavour, the only Vn and v T components at 
production are secondary neutrinos. However, the flavour mixing that takes places via the 
level crossings serves to even out the distribution across the three flavours. The final spectra 
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Figure 4. The final muon-flavour (anti)neutrino energy spectra, per annihilation, after propagation 
from the Sun's centre to Earth, for M x — 200, 500, 1000 GeV. The curves denote different annihilation 
branching ratios: equal branching to all lepton flavours (black, bold), branching only to electron 
flavour (blue, dashed), branching evenly to /z and r flavours (red, bold dotted). Also shown for 
comparison are the spectra arising from DM annihilation directly to t + t~ (green, dot-dashed) and 
to W+W- (brown, dotted). We set 9 13 ~ 8.8°, and (i = 1.2. 



have a similar shape, irrespective of the branching ratios chosen, with magnitudes that vary 
by only a factor of < 2. 

Scattering and absorption have a strong effect on the spectra for M x = 1000 GeV, but 
very little effect for M x = 200 GeV. This is expected as absorption depends exponentially 
on the neutrino energy. For large M x , the high energy regions of the spectra are strongly 
damped, removing the high energy peak. This is a bigger effect for neutrinos than for 
antineutrinos due to the different values of £ — the scattering in the Sun affects neutrinos 
more strongly than their antiparticles. Overall, absorption significantly reduces the number 
of detectable neutrinos produced per annihilation. For the energies we consider, regeneration, 
which reintroduces neutrinos creating a pile-up at lower energies, is a smaller effect. 

These plots in figure 4 were produced assuming a value for #13 of 8.8°. For the purposes 
of comparison, we also computed spectra for #13 = 0.1°. In the #13 ~ limit, the evolution 
through the high energy #i3-driven resonance is completely non-adiabatic, with Pq ~ 1 for 
all energies of interest, while for the larger mixing angle the resonance is (partially) adiabatic 
for energies < 1000 GeV. In figure 5, we compare the neutrino spectra for the two choices of 
^13. Decreasing #13 decreases the adiabaticity of the resonance, and increases the differences 
between the spectra for each flavour. This effect is most significant at low energies, where the 
level-crossing probabilities are small (see figure 3), while at high energies Pq — > 1 irrespective 
of the value of #13. For antineutrinos, there is very little sensitivity to #13 as the ^13-resonance 
condition is never met. From here on, we shall assume equal branching to all three flavours 
(/e = fn = /r)j which renders the results insensitive to the value of #13. 
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Figure 5. The spectra for M x = 500 GeV, for different choices of of 6*13, and where \x = 1.2. 
For branching only to e-flavoured leptons the spectra are shown for 6*13 = 0.1° (purple, solid) and 
9 13 = 8.8° (purple, dotted). For branching to fi and r flavours equally the curves correspond to 
9 13 = 0.1° (green, dashed) and 9 13 = 8.8° (green, dot-dashed). Note that there is no significant 
dependence on 6*13 for the case of equal branching to all flavours. 



4 Signals at Earth 

We now compare the signal from EW bremsstrahlung annihilation, to that arising in models in 
which DM annihilates instead to 2-body final states such as W + W~ or t + t~. These channels 
have been chosen for comparison as they both have relatively hard neutrino spectra, with a 
significant flux at large values of x. As a result the existing limits on these channels [16-18] 
are more stringent than those with softer neutrino spectra. As seen in figure 4, the EW 
bremsstrahlung, W + W~ and r + r _ channels all have broad energy spectra, of similar shape 
and scale. However the EW bremsstrahlung channel has a significant high energy peak at 
x ~ 1, which will improve its detectability compared to the W + W~ mode for which limits 
are commonly determined. We estimate relative event rates in the following subsection. 

The spectra in figure 4 have been computed per annihlation; the absolute normalisation 
of the flux is specified by eq. 2.1. (If we assume capture-annihilation equilibrium, the abso- 
lute normalization of the event rate is set by the capture rate alone times the annihilation 
branching ratio of the channel of interest.) However, for a given annihilation rate (which we 
leave as a free parameter) it is useful to compare the ratio of event rates for possible annihi- 
lation channels. We can thus assess the detectability of the EW bremsstrahlung annihilation 
channel relative to that for the comparison 2-body final state channels, W + W~ and t + t~, 
by comparing their relative event rates. 

4.1 Muon event rates 

We consider muon event rates in a detector such as IceCube. Although IceCube is able 
to detect neutrinos of all flavours via cascade events, we consider here only the muon-track 
events as they have the greatest detection probability, since the long muon range in ice allows 
the detection of muons created well outside the detector. 

We estimate event rates following a similar procedure to that in refs. [55, 67-69]. A 
muon neutrino can undergo charged current scattering with a nucleon within the ice of the 
detector, to produce a muon with energy = E u (l — y), where y is the charged current 
inelasticity parameter. This muon will then propagate through the ice and lose energy, with 
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a range given by 

" a + f3E u 



R ll {E„E^) = - Q \n 



1 



a + f3E^ 



(4.1) 



where a = 2.0 Me V cm /g, (3 = 4.2 x 10 _6 cm 2 /g, and E^ is the muon energy detection 
threshold, taken to be 100 GeV. Therefore, the probability that a neutrino with energy E u 
creates a muon that is detected with an energy above E^ r is given by 

P(E V , Ef r ) = p N N A a nucleon {E v )R^E v {l - y), E^ T ) , (4.2) 

where pw is the target nucleon density, N A is Avogadro's number, and cr nuc \ eon (E u ) is the 
charged current scattering cross section in ice. Finally, the muon track event rate is propor- 
tional to 

d -^P{E^E^)A^dE^ (4.3) 

where A eS is the effective area of the detector; for IceCube this is ~ 1km 2 [70]. The scattering 
cross section cr nuc i eon (E u ) is proportional to E u and, in ice, the neutrino cross section is 
approximately twice as large as that for antineutrinos. We choose the average values of 
y as tabulated in ref. [67], which are approximately ~0.45 and ~0.35 for neutrinos and 
antineutrinos respectively. As a result, the lower bound on the (anti)neutrino energies that 
will produce observable muons is 

Eu. > <7(1 - Ifr), Ep„ > <V(i - v„), (4 4) 

~ 182GeV, ~ 154GeV . 

In table 4.1, we show the ratio of muon-track events for various annihilation modes, 
assuming the same overall normalisation. For example, by setting the branching ratio in the 
table equal to 1, we can compare the event rate for a model in which the annihilation proceeds 
via bremsstrahlung with that for a model in which the annihilation is solely to W + W~ or 
t + t~ , assuming the same total annihilation rate. 

We see that, for a given annihilation rate, annihilation via the EW bremsstrahlung 
channel would produce more events than annihilation via a W + W~ channel, for all masses 
considered, with a ratio of ~ 4 at M x = 200 GeV. The ratio is largest for low dark matter 
masses, because only the high energy regions of the spectra are above the detection threshold, 
and hence the high energy peak in the EW bremsstrahlung spectrum dominates. The ratio 
is smaller for large M x , as the detection threshold then corresponds to smaller values of x, 
and absorption effects damp the spectra at high energy. There is only a mild dependence of 
these results on the value of the mass degeneracy, p = M^/M^. For larger values of p the 
high energy peak is diminished, but more mid-to-high energy neutrinos are produced, such 
that the total number of detectable events remains approximately constant. 

Detection prospects for annihilation channels which produce hard spectra, such as 
W + W~ , are much better than those for annihilation channels which produce much softer 
neutrino spectra, such as qq. Given that annihilation via the bremsstrahlung channel would 
result in more events than annihilation to W + W~ , its detection prospects appear promising. 
The experimental sensitivity for a given annihilation channel is usually expressed in terms 
of a limit on the DM-nucleon scattering cross section which controls the absolute event rate 
(assuming capture- annihilation equilibrium). For annihilation to purely VF + VF~, the pro- 
jected IceCube sensitivity (180 days of the full 86 string configuration) to spin-dependent 
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H = 1.2 


leV 
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500 C 
/x = 1.2 


fj, = 5 


1000 ( 
/x = 1.2 


3eV 
fj, = 5 


<3>„(brem) / Br(brem) 
<&v(t+t-)I Br(r+T-) 


1.8 


2.1 


1.2 


1.2 


0.6 


0.7 


<E>„(brem) / Br(brom) 
$ V {W+W-)I Br{W+W~) 


3.8 


4.6 


2.4 


2.6 


1.9 


2.0 



Table 1. The ratio of fi muon-track events in an IceCube-like detector, for various annihilation 
channels. We have assumed the same overall normalization for each possible annihilation mode, that 
is, assuming capture-annihilation equilibrium to hold, we assume the same capture rate for each case. 
For the EW bremsstrahlung annihilation channel, equal branching ratios to all lepton flavours was 
assumed (/ e = / M = / T ). 



DM-nucelon scattering is oTjy ~ few x 10~ 5 pb [71] for the DM mass range we consider. 
Since annihilation purely via the bremmstrahlung channel will produce more events than the 
W + W~ channel by a factor of a few, we estimate a sensitivity of a^® ~ 1 x 10~ 5 pb for this 
channel (see also ref. [72]). 

4.2 Parameter constraints 

We have presented spectra per annihilation, in the spirit of ref. [9], choosing to leave the ab- 
solute normalization as a free parameter. However, the detection prospects for a given DM 
model depend on the absolute event rate, which is determined by eq. 2.1 and eq. 2.2 and, 
in general, will depend on both the DM-nucleon scattering and the DM annihilation cross 
sections. However, importantly, if capture-annihilation equilibrium is established, then the 
event rate is controlled by the DM-nucleon scattering cross section alone (modulo astrophysi- 
cal uncertainties in the local dark matter density and velocity). We now discuss whether this 
condition is expected to hold in the DM models of interest for electroweak bremsstrahlung. 

Expressions for the capture and annihilation rate are given, e.g., in ref. [46]. In order 
for equilibrium to be established, the product of (i) the DM-nucleon scattering and (ii) the 
DM annihilation cross section must be sufficiently large. If one assumes a spin dependent 
DM-nucleon scattering cross section of ~ 10 -5 pb (the IceCube sensitivity estimated above), 
an annihilation cross section of (av) > W~ 2 {av)th is required for capture-annihilation equi- 
librium to be reached within the lifetime of the solar system, where (av) th ^ 2 x 10~ 26 cm 3 /s 
is the annihilation cross section expected for a thermal relic [73]. We discuss each of these 
cross sections in turn. 

The bremsstrahlung annihilation cross section is expected to be approximately two or- 
ders of magnitude smaller than that for a thermal relic, (o"f)i :)rem ~ 10- 2 (av) th . 2 This 
is because the annihilation cross section at thermal freezeout is dominated not by the 
bremsstrahlung channels, but instead by annihilation to 2-body final states. The 3-body 
bremsstrahlung processes dominate the annihilation in the late universe, because the 2-body 

2 Larger annihilation cross sections would imply either a subdominant DM component, or that the DM was 
produced not by thermal freezeout but by some more exotic mechanism. A number of recent papers have 
considered indirect detection and set upper limits on the EW bremsstrahlung annihilation cross section based 
on observations of antiprotons and gamma rays. For instance refs. [30, 39] both find limits of approximately 
(av) < 1CP 24 — 10 _25 cm 3 /s over the mass range 100-300 GeV. Ref. [42] performs a search for a bremsstrahlung 
annihilation signal in the Fermi LAT gamma ray data and sets bounds of approximately (av) < 10 _27 cm 3 /s 
over the same mass range. 
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final states are p-wave suppressed by a factor v 2 ~ 10 -6 ; in the early universe v is much 
larger and the 2-body final state is not strongly suppressed. See the discussion in section 4 
of ref [42] for further details. 

For a purely leptophilic model, the DM couplings to quarks arise at loop level. For the 
coupling of eq. 2.3 the DM-nucleon interaction arises at 1-loop level through a coupling to 
the Z boson, resulting in an effective interaction ([48]) 

Q 

£a = j I 2{xi fl i5x){qi f ii5q), (4.5) 

where Q is the loop-induced form factor. The resulting spin-dependent DM-nucleon cross 
section is shown in figure 2(a) of ref. [48], and can take values < 10 -6 pb. 

Therefore, for the purely leptophilic model described by the coupling of eq. 2.3, the 
product of the DM-nucleon scattering and DM annihilation cross sections are sufficiently 
small that capture-annihilation equilibrium would not hold (although we would be close to 
this limit, with event rates suppressed by a factor of less than 10 for some parameters). 
However, if the model were not purely leptophilic, the DM-nucleon cross section could easily 
be larger, such that capture-annihilation equilibrium would hold. For instance, couplings 
identical to that we consider arise in SUSY models. Embedding the toy model of eq. 2.3 into 
a full theory (e.g. the MSSM, or other theory beyond the standard model) is a more realistic 
scenario. In that case it is likely that several annihilation channels would be open, which 
would need to be combined with model dependent branching ratios. In this paper we take a 
purely phenomenological approach to study the EW bremsstrahlung signal in isolation. 

5 Conclusions 

High-energy neutrinos produced by the annihilation of dark matter captured in the Sun offer 
a useful dark matter search strategy. By comparing the predicted neutrino fluxes with future 
detections, we can potentially distinguish between different DM annihilation channels, thus 
providing information to discriminate among competing DM models. 

DM annihilation is usually assumed to be dominated by annihilation to 2-body final 
states. However, there are important DM annihilation processes with 3-body final states, 
namely bremsstrahlung processes in which a W, Z or 7 is radiated. The electroweak 
bremsstrahlung processes XX ~^ f fZ and XX f y W have been shown to be the dominant 
DM annihilation mode in certain popular models in which the lowest order annihilation pro- 
cesses are helicity suppressed. In the context of solar WIMP annihlation, these electroweak 
bremsstrahlung annihilation channels produce large, energetic, fluxes of neutrinos. 

We have computed the neutrino energy spectra resulting from DM annihilation via the 
electroweak bremsstrahlung channel. This consists of primary neutrinos produced directly 
in the annihilations, which have a distinguishing peak near the endpoint E v = M x , and 
secondary neutrinos produced in decays of annihilation products. By comparing the event 
rates for the electroweak bremsstrahlung annihilation channel with those for annihilation 
instead to W + W~ , we have shown that the EW bremsstrahlung channel has the larger event 
rate per annihilation and hence promising detection prospects. 
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